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Abstract. We present the current status of hybrid approaches to describe heavy ion collisions and their future 
challenges and perspectives. First we present a hybrid model combining a Boltzmann transport model of hadronic 
degrees of freedom in the initial and final state with an optional hydrodynamic evolution during the dense and 
hot phase. Second, we present a recent extension of the hydrodynamical model to include fluctuations near the 
phase transition by coupling a chiral field to the hydrodynamic evolution. 



1 Introduction 

Heavy ion collisions are an excellent tool to study the prop- 
erties of dense and hot nuclear matter 1 1 1. One of the most 
interesting features to explore in heavy ion collisions is 
the phase diagram of quantum chromodynamics (QCD). 
At high temperatures and large baryonic densities the de- 
grees of freedom of strongly interacting matter are sup- 
posed to be partonic. There are indications that this state 
of matter, the quark-gluon plasma, was created in heavy 
ion collisions at RHIC and behaves as an almost ideal fluid 
f2 3 4 5 1 . It is one of the major challenges to learn about the 
quark-gluon plasma (QGP) from the final hadronic state 
observed in the detectors. The phase transition of QCD is 
twofold: besides the transition from deconfined partonic to 
confined hadronic matter, the chiral symmetry gets spon- 
taneously broken at lower temperatures while the symme- 
try is restored in the high temperature phase. Lattice QCD 
studies have revealed that these two aspects form one crossover 
transition at least for vanishing baryochemical potential at 
critical temperatures given between Tc - [145, 165] MeV 
M6I7I . Studies of effective low energy models of QCD show 
that the phase transition is a first order transition at high 
baryon densities and low temperatures |8|. Consequently 
the first order phase transition line terminates at a critical 
end point (CEP) in the T - fiB plane of the phase diagram. 

A broad variety of experimental signatures is proposed 
to study the QCD phase diagram in experiment. One very 
famous signature is the 'horn' structure in the excitation 
function for the kaon to pion ratio [9|. It was attributed 
to the onset of deconfinement at SPS energies ifTOll . Other 
proposed signatures cover aspects of the chiral phase tran- 
sition, a suggested critical point 1 11 12 1 and the first order 
phase transition 1131 . It is, however, not clear how much 
of a potential signal is developed in a realistic scenario 
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of a heavy ion collision. The fast dynamics present dur- 
ing the creation and expansion of the fireball of strongly 
interacting matter makes it difficult to defer results from 
thermodynamic calculations directly. A thorough theoreti- 
cal understanding of phase transitions in the environment 
of a heavy ion collision is necessary to make profound pre- 
dictions for the running low energy RHIC program and the 
upcoming FAIR project. 

We present results from a hybrid model approach to 
heavy ion collisions 1 14 15 16 17 1 and give an outlook on 
chiral fluid dynamic models which allow for the explicit 
propagation of critical fluctuations through the phase tran- 
sition. 



2 The hybrid approach 

The hybrid approach describes a heavy ion collision in 
three steps, the inital state, the hot and dense stage with the 
phase transition and the final interactions after hadroniza- 
tion. By fixing the inital state and the freeze-out of particles 
from the hydrodynamic ally propagated densities one can 
focus on the effects of the phase transition and viscosity in 
the intermediate state. The model presented here is an ex- 
tention of the standard transport UrQMD model by the op- 
tional use of an intermediate hydrodynamic evolution [14]. 
This hybrid approach is realized in the latest version 3.3 of 
UrQMD and can be downloaded from www.urqmd.org. 

2.1 The initial state 

The inital collisions are obtained from the Ultra-relativistic 
Quantum Molecular Dynamics model ||I8|I9i When the 
Lorentz-contracted nuclei have passed through each other 
all nucleons have potentially interacted at least once. This 
is the earliest time to assume local thermalization, which 
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is needed for the validity of a hydrodynamic description of 
the expansion of the created matter. The transition time to 
the hydrodynamic evolution is 

- 27? 

fstart — , - ; W 

Vrcm - 1 

where R is the radius of the nuclei and ycm is the Lorentz 
gamma factor of the two colliding nuclei in the center of 
mass frame. To map the particles from the UrQMD to the 
hydrodynamic fields they are described by Gaussian dis- 
tributions with a width of cr = Ifm. The net-baryon den- 
sity, the energy density and the initial velocity profiles are 
then transformed onto a 3-dimensional space-grid to ini- 
tialize the hydrodynamic evolution. Thereby the event-by- 
event fluctuations in the initial conditions are taken into 
account. The spectators are not considered in the hydrody- 
namic evolution but propagated further in the cascade. 

2.2 The equation of state in the hydrodynamic 
evolution 

The most important input to the hydrodynamic expansion 
of the system besides the initial conditions is the equation 
of state. It strongly influences the dynamics of the sys- 
tems and naturally describes the phase transition in local 
equilibrium. In the framework of the hybrid model differ- 
ent equations of state can be implemented and their in- 
fluence directly seen as the inital and final state remain 
unchanged. The hadron gas equation of state, for exam- 
ple, allows for a direct comparision of the underlying dy- 
namics of the hydrodynamic treatment versus the transport 
sitmulation of the hot and dense region, since it includes 
the same degrees of freedom as in UrQMD and does not 
exhibit a phase transition. It can be used for baseline cal- 
culations to explore the influence of vanishing viscosity 
and local thermalisation. Other available equations of state 
are the Bag Model and the chiral equation of state. The 
Bag Model equation of state has a strong first order phase 
transition for all baryochemical potentials with a large la- 
tent heat. The chiral equation of state is obtained from a 
chiral hadronic Lagrangian that includes all baryons from 
the lowest flavor-SU(3) octet and the multiplets of scalar, 
pseudo-scalar, vector and axial-vector mesons. Recently a 
deconfinement equation of state has been derived by in- 
cluding partonic degrees of freedom in the chiral Lagrangian 
II20I . very much in the same manner as the Polyakov-loop 
extention of quark meson models. It has been successfully 
implemented to include investigations on the confinement- 
deconfinement transition to the framework of the hybrid 
model |21 1. 

The hydrodynamic expansion is performed by a full 
(3 + l)-dimensional SHASTA code r22'23'| solving for en- 
ergy and momentum conservation of an ideal fluid. 

2.3 Freeze-out 

When the system has expanded and diluted to a point where 
the hydrodynamic al prescription becomes questionable and 



local equilibration can not be assumed any longer, a tran- 
sition to a Boltzmann transport description is performed. 
In order to transfer the hydrodynamic densities to particles 
the Cooper-Frye freeze-out equation 

= / (2) 

is employed. It connects the boosted phase-space distribu- 
tion f(x,p) and the momentum-space distribution dN/d^p 
along the hypersurface cr. For the present calculations a 
gradual transition is chosen, which corresponds approxi- 
mately to an iso-eigentime transition. Transverse slices of 
thickness Az = 0.2fm are frozen out when the energy den- 
sity £ in each cell of that slice has dropped below a critical 
energy density e^nt- It is important to note that for the cor- 
rect transition the degrees of freedom of both sides, the 
hydrodynamic and the UrQMD |24|, need to be the same. 
Therefore, the last step in the hydrodynamic evolution is 
performed with an equation of state of a hadronic gas. The 
Monte Carlo sampling implementation of the Cooper-Frye 
equation conserves energy,baryon number, electric charge 
and strangeness in each single event. 



3 Results from the hybrid model 

We present results from the hybrid approach to heavy ion 
collisions, particle multiplicities and ratios of strange par- 
ticles to pions. 

3.1 Particle multiplicities and mean transverse 
mass excitation function 

To study the dependence of the equation of state, viscous 
effects and nonequilibrium dynamics the hybrid model re- 
sults can be compared to pure UrQMD 2.3 calculations. 
In Fig. ([T]i the multiplicities (top) and the mean transverse 
mass excitation function (bottom) are shown for pions (left) 
and kaons (right) compared to data from most central Au-l-Au/ 
Pb+Pb collisions at intermediate center-of-mass energies 
= 2 - 18A GeV. While the yield of pions is reduced 
in the hybrid model calculation compared to the transport 
calculation due to entropy conservation in the hydrody- 
namic expansion, as seen in the upper left plot of Fig. ([T]i, 
it does only show a weak dependance on the equation of 
state. In the kaon multiplicities, upper right plot of Fig. 
(|T]l, one clearly sees that the UrQMD 2.3 underpredicts the 
kaon yields in the whole energy range. In the Cooper-Frye 
freeze-out prescription the strange particles are produced 
according to their thermal distribution which leads to an 
enhanced yield in the hybrid model. Except for the Bag 
Model equation of state one sees a clear excess over the 
pure transport calculation, and a splitting between differ- 
ent equations of state. Since the transverse mass spectra 
are (Fig. ([T] bottom) more sensitive to the pressure in the 
transverse plan it shows a more pronounced dependence on 
the equation of state. For both pions and kaons the chiral 
and the hadron gas equation of state give similar results. 
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2. Excitation function 
strange-baryon to pion 



Fig. 1. Pion (left) and kaon (right) multiplicities (top) and mean 
transverse mass excitation function for central (h < 3.4fm) 
Au+Au/Pb+Pb collisions. The experimental data (full black sym- 
bols) are compared to hybrdi model calculations with different 
equations of state (lines with symbols) and pure UrQMD 2.3 cal- 
culations (dashed line) I2I5I26I27I28I29I 



Only the Bag Model equation of state yields lower pion 
mean transverse masses at higher energies as it is expected 
for a first order phase transition. Again in the Bag Model 
equation of state the values of the mean transverse masses 
for kaons, lower right plot of Fig. ([T]i, are even below the 
nonequilibrium dynamics of the transport model. 



for kaon to pion ratios (left) 
ratios (right) for most central 
Au-nAu/Pb-l-Pb collisions in the energy range ^fs = 2 - 200A 
GeV. The hybrid model calculations are the solid and dashed 
black line for two different freeze-out criteria. UrQMD 2.3 results 
and experimental data I2l5l26l27l28f29l are shown for compari- 
sion. 

propagate fluctuations directly in the hydrodynamic evo- 
lution. In models of chiral hydrodynamics [30 3 1 1 critical 
fluctuations at the phase transition can be studied directly 
and in non-equilibrium. For this purpose a field theoretical 
model that exhibits a phase transition is coupled to a back- 
ground hydrodynamic evolution to simulate the expansion 
and cooling of the fireball in the course of a heavy ion 
collision. An effective model of QCD that exhibits a chi- 
ral phase transition and a critical point is the quark meson 
model with constituent quarks |8 1 



(3) 



3.2 The kaon to pion ratio 

The strangeness production is enhanced in the hybrid model 
due to the local thermal equilibrium asssumed for the hy- 
drodynamic expansion |25|. We can, therefore, expect to 
obtain different results for the kaon to pion and strange 
baryon to pion ratios. These are shown for central Au-i-Au/Pb-i-Pb 
collisions at center-of-mass energies ^[s - 2 - 200A GeV 
in Fig. (|2|. There are two different freeze-out criteria of 
fcrit — 4eo (solid line) and fciit = 5eo (dashed line). The 
gray line is a pure UrQMD 2.3 calculation that underesti- 
mates the data in all four cases for higher energies. This is 
due to an over-production of pion yields compared to the 
strange particle production. As expected the hybrid model 
gives a better description of the data for energies above 
~ 5 GeV. At lower energies the assumption of local 
thermal equilibrium might not be justified and nonequilib- 
rium propagation seems to be necessary. 



+ -{d,nd^n)- 



U(o-, n) , 



with the constituent quark field q - (u, d), the couphng g 
between the quarks and the chiral fields (cr, n). The poten- 
tial is given by 



A 



(4) 



where the vacuum expectation of the sigma field is (cr) = 
f„ - 93 MeV and of the pion field is {n:} - 0. Iiqcr = f^ml 
with m„ - 138 MeV is the explicit symmetry breaking 
term. Then - - nm^/A^. A^ = 20 yields a sigma mass 
= 2A?f^ + a 604 MeV. In order to have zero poten- 
tial energy in the ground state a term Uo - m*/(4A^)-f^ml 
is added. 



4 Critical fluctuations in chiral 
hydrodynamics 

While the hybrid model gives an excellent framework to 
study the effects of an (local) equilibrium versus a full 
nonequlibrium transport propagation, it is not possible to 



4.1 The field equation of motion 

The order parameter of chiral symmetry breaking is the 
sigma field. From thermodynamic considerations it is ex- 
pected to fluctuate largely at the critical point. Our purpose 
is, however, to propagate the chiral fields in nonequilib- 
rium and study the effects of the phase transition. It is ex- 
pected that critical fluctuations at the critical point are di- 
minished while nonequilibrium fluctuations at the first or- 
der phase transition are enhanced. We will in the following 
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concentrate on the evolution of the sigma field and neglect 
fluctuations of the pion fields around their vacuum expec- 
tation value. The sigma field interacts with the quarks and 
antiquaries, which we assume to be in local thermal equi- 
librium. Due to this interaction the sigma field is damped 
and according to Langevin dynamics gets random kicks 
from the heat bath. The classical equation of motion for 
the chiral fields is 



5„5^o- -I- — + 2g^d„cr 
60- ^ 



d^; 



1 



(Iny E 



np(E) + T]d,cr^^, (5) 



where np{p) is the Fermi-Dirac distribution, dg = 12 the 

degeneracy factor and E - -^Jp^ + g^cr^ the energy of the 
quarks. The mass of the quarks is dynamically generated 
by a nonzero value of the sigma field. For the damping 
coefficient we take the value 77 - 2.2/ fm 132 j . Then the 
noise term is, in Markovian approximation. 



(6) 
(7) 



For the hydrodynamic evolution of the quark fluid the 
equation of state is needed. The pressure can be obtained 
from the thermodynamic relations in mean-field approxi- 
mation 



Veff(tr,;r, T) 



-2d J 



I 



J2n)'- 



■ log(l -He ^) + U{o-,n) 



(8) 

The strength of the phase transition at yUg = can be tuned 
by changing the strength of the coupling g. It is a crossover 
for small coupling and a first order phase transition for 
larger couplings, we use g = 5.5. Here, the effective po- 
tential has two degenerate minima allowing for the study 
of nucleation and spinodal decomposition I13I33I . For an 
intermediate coupling of gr = 3.63 a critical point with one 
very flat minimum can be found. Energy density and pres- 
sure are related by the thermodynamic relations 

p((r,n,T)^-V,«(cr,T) + U(cr), (9) 
dpicr, T) 

e(a,n,T)^T ' \l - p(cT,T) . (10) 
oT 

Besides the limited degrees of freedom in this model, the 
main difference to the chiral equation of state used in the 
hybrid model calculations above is that the field value is 
not the minimum of the thermodynamic potential. It is ex- 
plicitely propagated and, thus, the energy density and the 
pressure explicitely depend on the actual value of the sigma 
field. In the numerical simulations below the equation of 
motion of the sigma field (|5]l is solved by a staggered leap- 
frog algorithm. Initially we assume the sigma field to be 
in equilibrium with the fluid in order to be sensitive to the 
effects of the phase transition only. 



4.2 The hdyrodynamic expansion 

The equations of relativistic hydrodynamics of energy and 
momentum conservation are are 







(11) 
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Fig. 3. Time evolution of the intensity of sigma fluctuations at a 
first order phase transition. 



where Tfluid is the energy-momentum tensor of an ideal 
fluid. The source term reflects the energy dissipation from 
the field into the fluid. By using the explicit form of the 
equation of motion (j5]) we derive for the source term 



"t^^ field 



-(-g{qq),T - r]d,cr + ^)av (12) 



The inclusion of the source term guarantees energy and 
momentum conservation of the entire system. The fluid dy- 
namic part is again propagated by a SHASTA code. The 
energy density is initiated as the equilibrium energy den- 
sity at r = 160MeV, ellipsoidal in x-y-plane and homoge- 
nous in z-direction, smoothed by a Woods Saxon type dis- 
tribution. 



5 Numerical results 

The crucial quantity to look at is the intensity of the sigma 
fluctuations. It is given by 



dMr 



1 



(2ny2a>k {2ny2a)it 



(wf|o-t|- 



(13) 



This quantity can be considered as the number of sigma 
particles produced from the excited sigma field, once the 
nonlinearities can be neglected in the equation of motion 
of the sigma field. In order to take into account the mass 
change of the sigma particle in a hot environment the av- 
erage temperature T-^y in the hot region is calculated. With 
the sigma mass 



the energy of the modes is 



(14) 



(15) 



The time evolution of ( [T3| ) is shown in Fig. ^ for a first 
order phase transition and in Fig. Q for a critical point. In 
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a first order phase transition the intensity of sigma fluctua- 
tions between 4 and 11 fm is much larger than at a critical 
point. This is due to nonequilibrium effects at the first or- 
der phase transition where the effective potential has to two 
degenerate minima. Since we initialized the sigma field in 
equilibrium the increase in the sigma fluctuations is due 
to the phase transition. In the beginning of the evolution 
the sigma field is in the minimum which is stable above 
the critical temperatures T > Tc- It becomes the unsta- 
ble minimum for temperatures below the critical temper- 
ature T < Tc- The baiTier, which separates the two min- 
ima, however, prevents the sigma field from relaxing to the 
new equilibrium value. The damping to forces it further 
to remain in the local but unstable minimum. This eff'ect 
is known as supercooling. The noise term turnes out to be 
too small to make nucleation an effective mechanism for 
relaxation. It is more likely that the field will remain in its 
unstable configuration almost until the barrier disappears 
and it can relax via spinodal decomposition ||33l . 



6 Summary 

We have presented results from hybrid models to describe 
heavy ion collisions. The fully integrated Boltzmann and 
(3 -I- l)-hydrodynamic model is well established to test the 
effects of a nonequilibrium versus a (local) equilibrium 
evolution and of viscosity. By switching between different 
equations of state and comparing to the original UrQMD 
2.3 transport model the effects of a the phase transition 
have been studied. It is interesting to mention that the as- 
sumption of local thermal equilibrium in the hybrid model 
gives a good description of the 'horn' structure in /n'^ 
ratios without invoking a phase transition. The implemen- 
tation of a chiral fluid dynamic model, which combines the 
explicit propagation of the order parameter of the chiral 
phase transition and a fluid dynamic expansion of a quark 
fluid, enables us to study the phase transition in heavy ion 
collisions beyond the equation of state. We have included 
damping and noise term in the nonequilibrium dynamics 
of the chiral field in order to study the enhancement of 



the fluctuations at the phase transition and the relaxation 
to equilibrium. We have found that the intensity of fluctu- 
ations of the sigma field is significantly larger at the first 
order phase transition than at the critical point where large 
fluctuations in equilibrium would be expected. 



Acknowledgements 

This work was supported by the Hessian LOEWE initia- 
tive through the Helmholtz International Center for FAIR 
(HIC for FAIR). The computational resources were pro- 
vided by the Frankfurt Center for Scientific Computing 
(CSC). The authors thank Igor Mishustin and Hannah Pe- 
terson for fruitful and interesting discussions. 



References 

1. S. Schereref flZ., Prog. Part. Nucl. Phys. 42 (1999) 279. 

2. J. Adams et al. [STAR Coflaboration], Nucl. Phys. A 
757, 102 (2005) 

3. B. B. Back et al, Nucl. Phys. A 757, 28 (2005) 

4. I. Arsene et al. [BRAHMS Collaboration], Nucl. Phys. 
A 757, 1 (2005) 

5. K. Adcox et al. [PHENIX Collaboration], Nucl. Phys. 
A 757, 184 (2005) 

6. S. Borsanyi, Z. Fodor, C. HoelbHng, S. D. Katz, 
S. Krieg, C. Ratti and K. K. Szabo [Wuppertal- 
Budapest Coflaboration], JHEP 1009 (2010) 073 

7. W. Soldner [for the HotQCD collaboration], PoS LAT- 
TICE2010 (2010)215 

8. O. Scavenius, A. Mocsy, I. N. Mishustin and 
D. H. Rischke, Phys. Rev. C 64 (2001) 045202 

9. C. Alt et al. [NA49 Collaboration], Phys. Rev. C 77 
(2008) 024903 

10. M. Gazdzicki, J. Phys. G 23 (1997) 1881 larXiv:nucl-| 
th/9706036|. 

11. M. A. Stephanov, K. Rajagopal and E. V. Shuryak, 
Phys. Rev. Lett. 81 (1998) 4816 

12. M. A. Stephanov, K. Rajagopal and E. V. Shuryak, 
Phys. Rev. D 60 (1999) 114028 

13. I. N. Mishustin, Phys. Rev. Lett. 82 (1999) 4779 

14. H. Petersen, J. Steinheimer, G. Burau, M. Bleicher and 
H. Stocker, Phys. Rev. C 78 (2008) 044901 

15. R. Andrade, F. Grassi, Y. Hama, T. Kodama and 
O. . J. Socolowski, Phys. Rev. Lett. 97 (2006) 202302 
I arXi v:nucl-th /0608067 1 . 

16. D. Teaney, J. Lauret and E. V. Shuryak, RHIC," 
arXiv:nucl-th/01 10037 

17. T. Hirano, Prog. Theor. Phys. Suppl. 168 (2007) 347 
|arXiv:0704.1699 [nucl-th]]. 

18. S. A. Bass etal. Prog. Part. Nucl. Phys. 41, 255 (1998) 

19. M. Bleicher et al, J. Phys. G 25, 1859 (1999) 

20. V. A. Dexheimer and S. Schramm, Phys. Rev. C 81 
(2010) 045201 

21. J. Steinheimer, V. Dexheimer, H. Petersen, M. Ble- 
icher, S. Schramm and H. Stoecker, Phys. Rev. C 81 
(2010)044913 



EPJ Web of Conferences 



22. D. H. Rischke, S. Bernard and J. A. Maruhn, Nucl. 
Phys. A595(1995) 346 

23. D. H. Rischke, Y. Pursun and J. A. Maruhn, Nucl. 
Phys. A 595 (1995) 383 [Erratum-ibid. A 596 (1996) 
717] 

24. L. V. Bravina et al, J. Phys. G 25 (1999) 351 
I arXiv:nucl-th/98 100361. 

25. H. Petersen, M. Mitrovski, T. Schuster and M. Ble- 
icher, Phys. Rev. C 80 (2009) 054910 

26. S. V. Afanasiev et al [The NA49 Collaboration], Phys. 
Rev. C 66, 054902 (2002) 

27. L. Ahle et al [E866 Collaboration and E917 Collabo- 
ration], Phys. Lett. B 476, 1 (2000) 

28. J. L. Klay et al [E-0895 Collaboration], Phys. Rev. C 
68 (2003) 054905 

29. M. K. Mitrovski et al [NA49 Collaboration], J. Phys. 
G 32 (2006) S43 

30. I. N. Mishustin and O. Scavenius, Phys. Rev. Lett. 83 
(1999) 3134 

31. K. Paech, H. Stoecker and A. Dumitru, Phys. Rev. C 
68 (2003) 044907 

32. T. S. Biro and C. Greiner, Phys. Rev. Lett. 79 (1997) 
3138 

33. P. Chomaz, M. Colonna and J. Randrup, Phys. Rept. 
389 (2004) 263. 



